Bare electron dispersion from photoemission experiments 
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Performing an in-depth analysis of the photoemission spectra along the nodal direction of the high 
temperature superconductor Bi-2212 we have developed a procedure to determine the underlying 
electronic structure and established a precise relation of the measured quantities to the real and 
imaginary parts of the self-energy of electronic excitations. The self-consistency of the procedure 
with respect to the Kramers-Kronig transformation allows us to draw conclusions on the applicability 
of the spectral function analysis and on the existence of well defined quasiparticles along the nodal 
direction even for the underdoped Bi-2212 in the pseudogap state. 
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With modern angle-resolved photoemission spec- 
troscopy (ARPES) [jj, |2| one gets a direct snapshot of 
the density of low energy electronic excited states in the 
momentum- energy space of 2D solids H H H @. ah 
the interactions of the electrons which are responsible 
for the unusual normal and superconducting properties 
of cuprates are encapsulated in such pictures, but are still 
hard to decipher. One way to take into account these 
interactions is to consider electronic excitations as quasi- 
particles which, compared to the non-interacting elec- 
trons, are characterized by an additional complex self- 
energy • Extraction of the self-energy from experiment 
is thus of great importance to check the validity of the 
quasiparticle concept and understand the nature of in- 
teractions involved, but appears to be problematic since 
the underlying band structure of the bare electrons is a 
priori unknown. 

One can evaluate the interaction parameters taking the 
bare band dispersion from band structure calculations Q , 
however this unavoidably increases the uncertainty of any 
conclusions on the strength and nature of the interactions 
involved. A direct determination of the bare band struc- 
ture from experiment would be much more attractive in 
this sense. Previously, the bare band dispersion has been 
assigned to the high binding ener gy part of the experi- 
mental dispersion S]. In Refs. we have discussed 
that the bare Fermi velocity estimated from the nodal 
ARPES spectra using the Kramers-Kronig (KK) trans- 
formation is in reasonable agreement with band struc- 
ture calculations 0, an( i with an analysis of the 
anisotropic plasmon dispersion , although it has been 
pointed out that in order to quantify interaction param- 
eters such as coupling strength Q or self-energy E3 a 
precise and reliable approach of bare band determination 
is needed. 

In this Letter we introduce an approach to directly ex- 
tract the bare band dispersion and the self-energy func- 
tions from ARPES spectra. We show that the approach is 
self-consistent within the highest experimental accuracy 



available today, and, applying the procedure to the spec- 
tra from the underdoped Bi-2212, we demonstrate the 
validity of the quasiparticle concept even in the under- 
doped regime and in the pseudogap state. The detailed 
description of the procedure as well as its application to 
other samples of different doping level is given in Ref.[H. 

We start with a brief overview of the basics of the 
nodal spectra analysis within the self-energy approach 
illustrated in Fig. ^ Measuring the photoemission in- 
tensity as a function of the kinetic energy and in-plane 
momentum of outgoing electrons, 7(i^,k), one obtains 
access to the spectral function of the one electron removal 
which is supposed to reflect the quasipartical properties 
of the remaining photohole: its effective mass and life 
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FIG. 1: Bare band dispersion (solid line) and renormalized 
dispersion (red points) on top of the spectral weight of inter- 
acting electrons. Though intended to be general, this sketch 
represents the nodal direction of an underdoped Bi-2212. 
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time. These properties can be expressed in terms of a 
quasiparticle self-energy £ = £'+i£", an analytical func- 
tion the real and imaginary parts of which are related by 
the KK transform (the dispersion relation) 01 • Neglect- 
ing the effects of the energy and momentum resolutions 
as well as the influence of matrix elements 0, one can 
take / oc A(uj, k) , where ui = Ek — hv + <f> is the energy of 
the electron in the crystal with respect to the Fermi level, 
hv is the excitation energy and <f> is the work function. 
In turn, the spectral function can be formulated in terms 
of the self-energy: 
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A(u,-k) = — 



£» 



7T (w - e(k) - £'M) 2 + £"(w) 2 



(1) 



where e(k) is the bare band dispersion. 

In case there is no interaction, i.e. electronic excitations 
can live forever, the spectral function is a delta function 
with the pole u> — e(k w ) = and, e.g. for the nodal di- 
rection, can be represented by the solid line in Fig. ^ 
When interactions are present, the self-energy leads to a 
shifting and broadening of the non-interacting spectral 
function. The resulting picture is essentially that which 
is measured in ARPES (the blurred region in Fig. H] illus- 
trates this). If one neglects the momentum dependence 
of the self-energy, then, from Eq. QJ, the momentum dis- 
tribution curves (MDC(fe) = A(k) UJ=const ) have maxima 
at k m {w) determined by uj — e(k m ) — £'(w) = for a given 
to. In other words, £'(w) = to — s(k m ), is that which is 
illustrated in Fig. ^ by the double headed arrow. In the 
region where the bare dispersion can be considered as lin- 
ear (e = upfe), assuming weak k-dependence of £" along 
a cut perpendicular to the Fermi surface the MDCs 
exhibit a Lorentzian lineshape || with the half width at 
half maximum W and £"(w) = —vpW(uj). 

Thus, the determination of both the real and imaginary 
parts of the self-energy requires the knowledge of the bare 
dispersion (or, in the vicinity to Ep, an "energy scale", 
e.g., Vp 9]). The KK transformation gives an additional 
equation which relates these functions: 



£'(w) = - PV 

7T 
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x — UJ 
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where PV denotes the Cauchy principal value. This 
opens the way to extract all desired quantities from the 
experiment. 

In the vicinity of the Fermi level, one can determine the 
coupling strength as A = — (dY,'(uj)/duj) UJ =o — v f/ v r ~ 
1, where vr = (dk m /duj)~^ is the renormalized Fermi 
velocity. Differentiating the KK relation, 



duj 7T 

for an even £"(w), 
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A = PV 
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Using the above definition of the D operator, v f 
A 



v^ 1 - DVF, or 1 



1/Z, where 
Z = 1 - v R T>W 



(5) 



is the coherence factor (0 < Z < 1). 

In case W(ui) decays to zero or saturates on the scale 
covered by experiment, as it is expected for the scattering 
by phonons [18|, the parameters vp , A or Z can be eas- 
ily determined from the experimental values of vr and 
DW. In cuprates, however, the MDC width W along 
the nodal direction does not decrease or even saturate in 
the whole experimentally accessible energy region (up to 
uj m = 0.5 eV) and one can make only a rough estimation 
expanding DW = Dg m W exp + D£= m W mod , where W exp is 
the experimentally determined function of to, and W mo d 
is a model function which depends on both the high en- 
ergy cut-off, \lo c \ > \io m \, above which S"(w) starts to 
decrease or saturate and a model for these high energy 
tails. For a simple estimation one can take W exp = aio 1 
and W rnod = au? m which gives T>^ m W exp = T)™ m W mod , 
demonstrating that the contribution of unknown high en- 
ergy tails can be essential. 

Since without the knowledge of W mo d(oj) it is not pos- 
sible to determine precisely the energy scale of the bare 
band or interactions from just the low binding energy 
part of ARPES spectra, in the following we examine a 
wider energy range of the ARPES data in order to solve 
the problem. In this energy range a deviation of the bare 
dispersion from a line should be taken into account. In 
Ref. |9| we have shown that the shape of the bare band 
can be determined from the shape of the Fermi surface 
through the tight binding (TB) fitting procedure but the 
energy scale must be determined independently. 

Along the nodal direction the TB band in the occu- 
pied part can be well approximated by a simple parabola 
e(k) — u>o(l — k 2 /k F ), where the energy scale parameter 
is the bottom of the bare band, ujq, or the bare Fermi 
velocity, vf — —2ojo/kp. Then, using this dispersion in 
ifTjl. one can express the complex self-energy function in 
terms of two experimental functions k m {ui) and W{to): 



£» 



V F 
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The KK transform £" — -> £' completes the system but 
adds new parameters which describe the high energy tails 
of £'». 

As an example, we use the experimental data from an 
underdoped Bi(Pb)-2212 (T c = 77 K) measured along 
the nodal direction at 130 K with 27 eV excitation en- 
ergy, at which only the bonding band is visible Il2| . The 
experimental details can be found elsewhere [lj, [l5j . 

Fig. |21 shows the results of the procedure we describe. 
The real part of the self-energy is determined in two ways: 
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FIG. 2: Real and imaginary parts of the self-energy extracted 
from the experiment with the described procedure. A com- 
plete coincidence between the corresponding parts of the self- 
energy calculated from the two different experimental func- 
tions, the MDC dispersion and MDC width, demonstrates 
the full self-consistency of the ARPES data treated within 
the self-energy approach. 



experimental error. These deviations can be expressed 
for these two curves as a nondecreasing difference in the 
slope ratios below and above the 70 meV kink. To visu- 
alize this we fit the functions in the range of \u>\ < 
170 meV to £[ ow (lj) = — Xu — AX(uj — Wfc){arctan[a;fc/<5] — 
arctan[(cj — cjfc)/<5]}/7r (see Fig- EJ 7 where ujk « —70 meV 
is a kink energy and S « 30 meV is a kink width |l9j . 
For 6 <C \oJk\, A and A^ = A — A A reflect the slopes below 
(\uj\ < \u>k\) and above the kink energy respectively, and 
it is the difference in X/Xh between EJjisp ano - ^kk which 
is irreducible by varying the parameters luq, w c and the 
high energy tails. 

The origin of this difference is related to resolution 
which mainly effects the low energy part of T,' kk (uj) 
through W{uo). Therefore, high energy slopes Xh and 
high energy tails which are resolution independent should 
coincide for T,' disp and S^k after the fitting proce- 
dure is performed. This appeared to be true and 
the final result is shown in Fig. [5] It is remarkable 
that the difference E^k — S^ isp completely coincides 
with R'(lu), the contribution from the overall resolu- 
tion R = 15 meV, calculated as R' = KK R", where 
i?"H = vAR 2 + [S"(w) - S"(0)] 2 - [E"(w) - E"(0)] is 
the contribution of the overall resolution to the scatter- 
ing rate. Such a complete coincidence substantiates that 
the self-energy constructed using Eqs. © and JJJ is self- 
consistent within the experimental accuracy currently 
available with ARPES. In order to check the correctness 
of the KK numerics, we also plot Y/^ k {uj) function which 
is obtained by back KK-transform of Tf KK {uS). 

For underdoped Bi(Pb)-2212, for which k F = 0.471 



(i) £^ isp , by Eq. © with vp as a parameter; (ii) by 
Eq. (JJJ and subsequent KK transform of symmetrized 
^width — * ^kk with additional parameters which de- 
scribe the tails, viz. a cut-off energy, u> c , at which 
the tail is attached and a model of its decay (see for 
details). If the self-energy we are trying to extract satis- 
fies the KK relation, S^ isp and should be identical 
within experimental uncertainty which sets the space of 
the parameters in question. 

In the experimental energy range of 350 meV we have 
found that Vp is restricted by such a procedure in con- 
trast to cu c and to the model for the high energy tails, for 
which only a relation can be obtained. In practice, we 
checked several models where T l 'J nod (uj) saturates to some 
maximum value, slowly decays or drops to zero above ui c , 
and for each model the fitting of E^k to S^ isp results in 
the same vf (or ojq) but different uj c (for details see Fig. 2 
in and related discussion). 

At this stage we have found that although the fitting 
procedure is working well, i.e. the standard deviation has 
a minimum as a function of (luqi w c), and the result- 
ing curves Ej isp and E^k almost coincide (not shown), 
there are still deviations between them which exceed the 
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FIG. 3: The bare band dispersion along the nodal direction 
of an underdoped Bi(Pb)-2212 (black solid line) on top of 
its spectral weight at 130 K measured by ARPES. MDC (or 
renormalized) dispersion shown by solid red line. 
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A 1 and vr — 2.04 ± 0.05 eVA, we have obtained for 
the bare band dispersion along the nodal direction vf = 
3.82 ±0.17 eVA (w = -0.90 ±0.04 eV). This bare band 
is shown in Fig. [31 on top of the electron spectral weight 
as measured by ARPES. The interaction parameters thus 
arc the following: A = 0.87 ± 0.12 and X h = 0.24 ± 0.03. 
As we discussed above, the energy cut-off, ui c , depends 
on the model used for the high energy tails of E"(u;) but 
we can state that \u c \ » |w |/2. 

The linear behavior of S'(w) over a wide en- 
ergy range \uj\ < \uik\ indicates, using the criterion 
lim^^o E"(w)/o; = 0, the existence of well defined quasi- 
particles in the pseudogap state: for the underdoped 
Bi(Pb)-2212 at 130 K the coherence factor Z = 0.54 ± 
0.03. The offset of E"(w) not only comes from finite 
resolution but also finite temperature and scattering on 
impurities, which are mostly energy independent fl5| and 
do not contribute to the slope of E'(u;) and, therefore, to 
the coherence factor. 

In we have noticed that the scattering rate at room 
temperature looks more linear for underdoped samples 
than for overdoped ones that is in favour of the marginal 
Fermi liquid model (MFL) |2fJ . It is important to stress 
that E'(u;) determined with better accuracy exhibits a 
linear behavior below and above the kink energy u>k (see 
Fig. [21 which is now difficult to reconcile with the MFL 
model: as far as a slope in E'(u;), according to Eq. ©, is 
mainly determined by the coefficient at {uj — uj x ) 2 term 
in the expansion of E"(a;) around lu x , the straight sec- 
tions on E'(w) imply the regions where E"(w) is precisely 
parabolic (exhibits constant curvature over some finite 
energy regions). 

Another point arises as a consequence of the tight cor- 
relation between E' and E" . Recently we have shown 0] 
that two different channels can be distinguished in the 
scattering rate: the doping independent Auger like de- 
cay which originates from the electron-electron Coulomb 
interaction and the doping dependent one which can be 
naturally assosiated with spin excitations. While such 
a decomposition of the scattering rate into two channels 
seems to be becoming commonly accepted |2lj . there is 
still a controversy about the origin of the doping depen- 
dent one. The present analysis shows that regardless of 
the nature of this channel, its doping and temperature de- 
pendence should appear in the doping and temperature 
dependence of E' and, consequently, of the renormalized 
dispersion, although it is clear that the variations in the 
latter should be marginal. This casts doubt on the ap- 
parent universality of the nodal Fermi velocity j^, viz. 
that under close scrutiny it will not hold. 

In conclusion, we demonstrate for the first time the 



full self-consistency of the data obtained using angle re- 
solved photoemission and treated within the self-energy 
approach. This self-consistency is demonstrated on an 
underdoped Bi-2212 sample (T c = 77 K) in the pseu- 
dogap state. The extracted bare band dispersion is in 
good agreement with the band structure calculations 
and allows one to quantify the self-energy of the elec- 
tronic excitations in the real energy scale. The accu- 
rately determined real and imaginary parts of the self- 
energy prove the existence of well defined quasiparticles 
along the nodal direction even in the pseudogap state of 
Bi-2212. 
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